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SUMMARY 


The problem of predicting vortex breakdown, which is characterized by a drastic 
change in core size and trajectory, is still basically unsolved. In the present 
paper, viscous flow m a long rotating pipe is analyzed, as a model for the type of 
breakdown observed in wake vortices. The full Navier-Stokes equations are shown to 
simplify to a set of unlinked partial differential equations. These equations allow 
for two separate solutions, for both full pipe flows (the classic rotating Hagen- 
Poiseuille flow) and cavity flows, in which a central area, with no net mass flux or 
viscous dissipation, is produced. Conditions permitting switching from one type of 
flow to the other, with no change in the mass flux or velocities at the pipe sur- 
face, are found. It is then argued that these conditions indicate when vortex 
breakdown will occur. 


1 . INTRODUCTION 


The bathtub vortex, with its hollow center, was one of the first observed and 
best known flow patterns. Swirling flows of liquids m pipes can also produce 
configurations with hollow centers under certain conditions of mass flux and tangen- 
tial velocity (refs. 1,2). The free surfaces thus produced can support standing 
waves (ref. 2), and the flow field obtained is similar in some respects to that 
obtained when axisymmetnc vortex breakdown occurs m tubes. This similarity has 
already been noticed (refs. 3,4) and related to the possibility of obtaining an 
example of conjugate flow states. The existence of conjugate flows forms the basis 
of Brooke-Benjamin's theory of vortex breakdown (ref. 3). 

The studies cited above all deal with inviscid flow. In the present work, 
viscous flow will be examined to investigate the effects of viscosity on the produc- 
tion of hollow cores, as well as to better describe this form of pipe flow as a 
model for a free vortex core. Thus, we look at the viscous flow inside a rotating 
tube whose length is much greater than its diamater. This can be seen as an 
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idealization of a vortex core, if the coordinate system in which the tube wall is at 
rest is considered to be moving at the free-stream velocity. 

In section 2, the equations of motion are solved for the flow configuration 
above. Solutions showing cases of cavity flow arising for high swirl ratios are 
shown (sec. 3) to be a possible model for the inception of vortex breakdown of the 
axisymmetric kind m vortices of uniform strength. This approach is similar to the 
quasi-cylindrical, boundary-layer-type analysis (see ref. 5 for a review), in that 
conditions causing failure of the cylindrical approximation will indicate actual 
occurrences, such as separation or breakdown, even though the assumptions on which 
the analysis is based are thereby invalidated. 


2. ANALYSIS 


To obtain conditions for a change in the flow configuration, we now look at 
flows with fixed mass flux Q in a tube of fixed radius a. Considering a viscous 
incompressible fluid, the equations of motion for steady flow in fixed cylindrical 
coordinates (r, 9, x; see fig. 1) are, after nondimensionalizing , 
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Here p and v are the density and kinematic viscosity, respectively, and u, v, and 
w are, respectively, the axial, radial, and tangential components of velocity 
divided by the average speed 0 defined below. All lengths are normalized by the 
radius a. The pressure is divided by a reference gauge pressure p , at the 
entrance; 0 is defined by 



pita 


(2) 


The tube (which can describe the outer boundary of a vortex core) is rotating 
around its longitudinal axis at speed W. This can be normalized as 
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( 3 ) 
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The boundary conditions for equation (1) are, at r = 1 (the tube wall), 

u = 0 
v = 0 
w = W 

and at r = b (the cavity) where 0 < b < 1, 



v = 0 


(4a) 

(4b) 

(4c) 
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(5b) 


3w 

3r 


0 


(5c) 


This is reminiscent of annular Poiseuille flow with the inner cylinder behaving like 
a free surface. The area covered by r < b will be treated later. 

The classic case of rotating Hagen-Poiseuille flow is obtained for b = 0 in 
equation (5). The boundary conditions do not depend on x, so one can assume 

u t u(x) (6) 

Also, the boundary conditions exhibit axisymmetry so that 

38 = ° (7) 


in equation (1) for these cases. From equations (6), (7), the equation of continu- 
ity and boundary conditions (4b) and (5b) we obtain: 

v = 0 ; b < r < 1 (8) 
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and equations (1) reduce to 
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where R is the Reynolds number (R = Ua/v). 

From equation (6) we deduce that 3p/3x is a constant (which we will designate 
as P x ). Thus, equation (9a) can be solved independently of equations (9b) 
and (9c). The full solution for (9a) is 


2 

u = c^ Hn r + c 2 r + c^r + c Q 
Substituting equations (4a) and (5a) into (10) we obtain: 


( 10 ) 
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where 


A 



(12) 


is a pressure coefficient. 

One can obtain the nondimensional mass flux by integrating the velocity over 
the tube. This is now done for three cases. 

Case 1 . The classic Hagen-Poiseuille case (H-P flow) . This case is obtained 
when b = 0, and yields 


= r ( - p x )(1 - r2) 


(13) 


The nondimensional mass flux is then 

f 
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Case 2. No mass-flux in r<b(0<b<1). This is the case either when the 
boundary at r = b is a solid cylinder moving axially at a speed so that 
3u/3r =0 at r = b, or when a free surface is produced and maintained by rotation 
of the outer tube. It is the latter case that is of interest. Here, if the pres- 
sure is exactly the vapor pressure of the liquid, the internal hollow core is filled 
by vapor with much lower density, so that the mass flux in r < b may again be 
neglected. This can also describe cases of axisymmetric vortex breakdown with 
elongated closed "bubbles" being produced (ref. 6), for the net mass flux in a 
closed bubble is zero. 

The mass flux is now 


r 


Q„ = 2u J ur dr 
d b 


(15) 


with u taken from equation (11). Integration yields 


Q 2 = RA(-P x )(1 + 3b 4 - 4b 2 - 4b 4 an b) 


( 16 ) 


The third term is alwayys smaller than unity for 0 < b < 1 , and goes to zero when 

b = 1. Thus, for equal mass flux, the pressure gradient in this case has to be 

larger than that for the H-P flow in the same tube. If we define the pressure 
gradient ratio for a given fixed mass flux to be 

(-P ) H-P 
a ~ (-P x ) annulus 

we can obtain the following relation (shown in fig. 2) between this ratio and the 
cavity radius b: 

a = 1 + 3b 4 - 4b 2 - 4b 2 in b ( 18) 

Case 3. Uniform flow in r < b (internal plug flow) . This case is a good 

approximation to the situation in which part of the flow is turbulent (at r < b). 

fts such, the equivalent viscosity is much higher (eddy viscosity) in this part of 
the tube, and the velocity gradients much smaller. We thus approximate the flow 
in r < b to have a uniform axial velocity, obtained by substituting r = b m 
equation (11) 



(19) 


U b = IT ( -P x ) [2b 2 an b + 1 - b 2 ] 


The nondimensional mass flux is now obtained by 
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( 20 ) 


where u is taken from equation (11). After integration, we obtain 


Q 3 = RA(-P x )[1 - b 2 ] 2 


Again, defining the ratio of pressure gradients as s, 
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The variation of 6 with b also appears in figure 2. 

Returning to equation (9), we now solve equation (9c). This can be recast as 
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which is an Euler equation, with solutions 


w = — + c„r 
r 2 


(23) 


r = 


Applying boundary condition (4c), and requiring that the free surface at 
b be stress-free we have 


r = 1 ; w = W 


(24a) 


r = b ; |H = 0 

’ 3r 


(24b) 


as the boundary conditions that equation (23) has to fulfill. The result is 
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which gives, when b = 0, 


w = Wr 


as expected. Equation (25) can now be substituted into equation (9b) to obtain the 
radial pressure gradient: 
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The pressure difference in the radial direction is 


Z* 1 3D . A W 2 
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When b = 0 (no cavity) the pressure difference is 


Ap = £ W 2 


so that for 0 < b < 0.528 the pressure difference in the radial direction is 
greater in the case of cavity flow than for the regular rotating pipe flow, and 
for 0.528 < b < 1 it is less (see fig. 3). The nondimens ional stress on the outer, 
pipe wall m the tangential direction is * 


J_ aw 
T w ‘ R 3r 


W 1 - b 
1 1 + b 2 


where t is the shear stress acting on this pipe; therefore, the torque on the 
tube is less for all cavity sizes than for full pipe flow. 

The stability of the flow field described by equations (11) and (25) was 
checked, by means of the so-called Rayleigh, Synge, and Howard-Gupta criteria 
(ref. 7). No axisymmetric instability was found, so that cavity flow behaves like 
Hagen-Poiseuille flow in this respect. 
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3. COMPARISON WITH EXISTING WORK 


In the previous section, solutions were developed for viscous pipe flow with 
swirl allowing for an internal axisymmetric free surface. Increasing the axial 
pressure gradient for given tube mass flux and tangential velocity can cause the 
flow to develop an internal axisymmetric cavity, or cause an existing cavity to grow 
radially. This is in qualitative agreement with experimental results by Escudier 
et al. (ref. 1) who caused growth of the bubble diameter by putting a sudden con- 
traction in the pipe diameter (which causes an axial acceleration corresponding to 
an increase in the pressure gradient). Quantitative agreement is not to be expected 
here, for the experimental setup was somewhat different. In the experiment, there 
was rotating flow in a stationary pipe, and the tube m the theoretical model repre- 
sents the edge of the vortex core in the experiment. Nevertheless, as we shall see, 
various experimental observations can be reproduced in an encouraging manner. 

The radial pressure difference across the pipe is given by equation (27) and is 
shown in figure 3. As we can see, in the range 0 < b < 0.528 there are, for a 
given fixed mass flux, always two cavity diameters that result in the same radial 
pressure difference. As a specific example, take the case b = 0. This describes 
rotating Hagen-Poiseuille flow. A change in the axial pressure gradient, by a 
factor of 3.16 (fig. 2) results in a flow with a nondimensional cavity radius of 
b = 0.528, having the same radial pressure difference between the center of rotation 
and the tube wall. This means that removing the outside tube wall would not result 
in any change in the radial direction (so that this could model a viscous vortex 
core in a fluid instead of in a solid tube). 

The axial pressure drop along the pipe is the viscous analog to Brooke- 
Benjamin's "flow-force" (ref. 3). This is more aptly seen here as a kind of "avail- 
ability" in the thermodynamic sense because of the dissipation resulting from vis- 
cosity — as manifested by the longitudinal pressure gradient. Thus, the cavity 
(b = 0.528) flow has to lose more energy to traverse a given length. Having less 
available energy, the cavity flow can be seen as the subcritical conjugate in 
Brooke-Benjamin's terminology. Although the results above are mathematically rigor- 
ous, at the present time our view is that the analogy should be treated as an 
explanatory device only. 


4. ORIGIN OF CAVITY FLOW AND ITS RELATION TO VORTEX BREAKDOWN 


As mentioned m the Introduction, the pipe flow of case 2 (sec. 2) can serve as 
a model for flow in the vortex core, if the outside wall is considered to move at 
the speed of the potential flow external to the core. The internal mass flux is 
then the axial flow obtained in the vortex core (ref. 8). 
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For flow m the tube, the pressure at the outside wall diminishes in the down- 
stream direction along the tube axis (see fig. 4). It follows that at some point 
x c the pressure at the wall will be exactly equal to 

P c = |w 2 (31) 


When this point is reached, the pressure at the tube center vanishes (see eq. (29)); 
that is, a cavity can potentially be formed. Equation (31) can be used to establish 
a position at which the cavity solution is first obtained, if the pressure at the 
outlet is taken as zero. This position is a function of the inlet pressure (the 
reference pressure p m used in the definition of A), the mass flux, and the tan- 
gential velocity ratio. Written in terms of the present nondimensional quantities, 
this is 


<- p x )X c = P. - P 0 (32) 

The ambient pressure at the pipe outlet was taken to be zero so one can write 
p m = (-P x )L , where L Q is the length of a fictitious tube along which the total 
pressure excess is lost to friction. The length L Q is greater than any real tube 
that sustains a given mass flux. Substituting this relation, as well as equa- 
tion (31), mto (32) we obtain 


p 0 = <-P x )(L 0 - * 0 ) (33) 

For full pipe flow (rotating Hagen-Poiseuille flow) the pressure gradient (-P x ) can 
be obtained from equation (14). As a result, equation (33) can be written as 


2 2 

r _ RA^W 

o - x c 64 


(34) 


Length L Q has to be, as mentioned above, greater than the real tube length. 
Putting the actual tube length as x Q m equation (34) will thus define minimum 
values of the parameters for a cavity to exist. The point at which a cavity can 
first be sustained is thus a function of the swirl ratio W, the Reynolds number, 
and the mass flux. 

The next step is to attempt to relate the findings above to vortex breakdown, 
which is essentially a two- or three-dimensional phenomenon, whereas up to now the 
analysis was one-dimensional. As a result, we cannot expect predictions of the 
actual postbreakdown flow, but we suggest that the axial position of the breakdown 
is related to the critical point obtained above. Once zero pressure is obtained, 
any further drop m pressure is impossible. The constant axial pressure gradient 
required to fulfill equation (6) no longer exists, and the one-dimensional flow 
field has to disappear beyond this point. One can then use the boundary-layer 
separation analogy as proposed by Hall (ref. 5). It states: "If, at some location 
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m a quasi-cylindrical vortex core the results show appreciable instead of small 
axial gradients, there must also be appreciable axial gradients at that location in 
the corresponding real vortex core, even though the quasi-cylindrical approximation 
must fail there." On this view, x c in equations (33) and (34) may be seen as an 
approximation to the breakdown position, with an accuracy of 0(1/x c ). Our results 
can therefore only be qualitatively compared with those of Escudier and Zehnder 
(ref. 9). They forced breakdown at the entrance to the contraction in their 
apparatus, which is equivalent to the entry zone in the present model, that is, 
x c = 0. Escudier and Zehnder showed that for a fixed breakdown position, the 
Reynolds number, mass flux, and circulation can be related by means of an expression 
of the form (in present notation) 


RW y Q 6 = const (35) 

where y and 6 are empirically obtained values. As mentioned above, no quantita- 
tive comparison is possible (also a result of Escudier and Zehnder' s technique for 
producing swirl), but equations (33) and (34) nonetheless show that relationships 
such as equation (36) can also be obtained for pipe flow. 


5. CONCLUDING REMARKS 


Where it exists, the free cavity has been shown both experimentally and by an 
mviscid analysis to be capable of sustaining standing waves (ref. 2). This fact, 
in addition to the higher dissipation (resulting in lower available energy), indi- 
cates that the cavity flow obtained for a constant mass flux may be a viscous sub- 
critical flow m the sense of Brooke-Benjamin (ref. 3). Criticality conditions such 

as Squire's (see ref. 5) can be retrieved by assuming inviscid flow, but they are 

not really relevant here. 

We thus conclude that equation (35) can be used as an approximate predictive 
tool for breakdown m columnar vortices. Increasing vortex strength or swirl ratio 
causes x Q to dimmish, resulting m earlier breakdown. This situation can be 
qualitatively related, in the following sense, to that of vortices in the wake of 
slender wings or bodies. Increasing the angle of attack for a given configuration 
causes both vortex strength and swirl ratio to go up, with a corresponding forward 
movement of the point of breakdown. Increasing sweep angle at a fixed angle of 
attack causes W to decrease and R and L Q to increase. The point of breakdown is 
predicted to move rearward, and this is actually observed (ref. 10). 
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Figure 2.- Ratio of pressure gradients required to drive full, to annular rotating 
viscous pipe flow a and full laminar, to part turbulent flows S vs. ratio of 
annulus diameter to pipe diameter b. 



Figure 3.- Pressure increment in the radial direction, normalized by the pressure 
increment for full rotating pipe flow vs. ratio of annulus diameter to pipe 
diameter b. 
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Figure 4.- Schematic representation of pressure distribution and cavity location, 
(a) Distance to first cavity position x c , tube length L, and hypothetical 
maximum tube length L Q (eq. (33)) along which the total pressure excess is lost 
to friction, (b) Schematic pressure distribution. All lengths are normalized 
by a, the tube radius. 
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